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COORDINATED SCIENCE LABORATORY 
SUMMARY OF 

PROGRESS REPORT FOR JUNE, JULY, AUGUST, 1964 


1. _ Aerospace Group 

Three rockets employing closed loop instrumentation for 
radio propagation measurements in the D region of the ionosphere were 
fired before, during, and after sunrise on July 15, 1964, Reduced 
data from the April 16 shot is presented. 

2. _ Surface and Atomic Physics 

Results of tests of an ultrahigh vacuum, rotary-motion, 
feedthrough are reported. An experiment to study the elastic and 
inelastic scattering of electrons at solid surfaces is described. A 
study of the cross-section for scattering and ionization of barium 
atoms by low energy electrons is reported. A summary of the progress 
on experiments to study the adsorption-desorption kinetics of gases at 
solid surfaces, the angular distribution of secondary electrons, and 
Auger electron ejection is given. 

3. _ Computer Application Group 

Continuing developments on an experimental time-sharing 
system are reported. Improvements and additions to the SMP-CSX-1 
Bubble Chamber Data Processing work are described. 
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4. _ Systems 

An important problem in the computer analysis of networks, 
namely the generation of trees without duplication, has been solved. 

The necessary and sufficient conditions for the realizability of a 
given state diagram as a quasi-linear sequential machine have been 
obtained. The necessary and sufficient conditions for the realizability 
of a given matrix as the A matrix (state matrix) of an RLC half-degenerate 
network have been found. The extension of the Dasher S 3 mthesis pro¬ 
cedure into the right half plane is complete. Further progress is 
reported on the problems of the path matrix of a contact network, lossy 
communication networks, and self-diagnosing computers. New results in 
coding theory relate to the transmission rate in a communication channel. 
Several studies of nonlinear oscillating systems have been completed 
and are reported. In the area of control, progress is reported in 
studies of stability, system optimization, and rendezvous strategies. 

^ _ PLATO 

Progress on the construction and circuitry of the twenty 
student station teaching system indicates that at least ten student 
stations will be operable in November. The compiler, CATO (Compiler 
for Automatic Teaching Operation), has been completed. The resident 
program for CATO is now completely revised to allow for more efficient 
use of the PLATO III equipment, and compatibility with CATO-compiled 
programs. The inquiry training lesson, REPLAB, has been used further, 
this time by a demonstration class in a teacher's summer workshop on 
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inquiry training techniques. A short experiment in teaching small 
children the letters of the alphabet has suggested some of the problems 
involved in using PLATO with very young children. New lesson material 
under development includes a text book tester. The feasibility of 
externally controlled discharges has been successfully demonstrated in 
connection with the plasma discharge display tube research project. 

6. _ Vacuum Instrumentation 

Work on ultrahigh vacuum mass spectrometer studies is reported. 
Modifications to the photocurrent suppressor gauge show a considerable 
improvement in performance is possible. Some methods of repairing poor 
quality or damaged storage tubes are discussed. Impact ionization 
coefficients for CO on molybdenum and O 2 and CO on platinum are reported. 

7 . _ Plasma Physics 

An algebraic theory of shock structure has been completed. 

The results of this theory and of recent Monte Carlo calculations have 
led to the preliminary design of a new numerical integration scheme 
for the Boltzmann equation. Theoretical work on oscillations in non- 
uniform plasmas was started. One particular model is discussed in 
some detail. The eigenvalues of oscillation are calculated explicitly. 

8. _ High Magnetic Field Superconductors 

Progress is reported in the following areas; (1) Preparation 
of single crystalline Nb^Sn; (2) Studies of current peaks in the current- 
voltage characteristics of superconducting tunneling junctions related 
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to the AC Josephson effect; (3) Superconducting properties of thallium- 
antimony alloys; and (4) Development of a very sensitive low impedance 
voltage detector. 

9. _ High Voltage Breakdown 

The significance of critical field values for breakdown is 
discussed. Observations of emission points on broad electrodes, the 
effect of gas on prebreakdown currents, the use of gold electrodes, and 
further observations of field-emission flicker are reported. 

10. _ Thin Films 

Photoconductive measurements have been made on tin oxide 
films in the continuing effort to learn more about the electronic 
structure of this material. Using heated substrates, films of vanadium 
and niobium have been made which have normal superconducting properties. 
This is an important improvement over earlier results. The substrates 
were heated using an evaporated titanium film on the reverse side on 
the substrate as a resistive heating element. Electron diffraction 
and electron microscope photographs have been made of the structure of 
some tin films grown on cleaved NaCl substrates. Preliminary attempts 
have been made to deposit by evaporation luminescent screens. 

11. _ Computer Operations 

Operating statistics, equipment modifications, and systems 


programming are reported. 
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1, AEROSPACE GROUP 

H. Knoebel H. Krone W. Prothe 

D. Skaperdas J. Gooch R. Anderson 

B. Kirkwood G. Karr 

On July 15, 1964, the Coordinated Science Laboratory fired 
three rockets of the Nike-Apache type in a series of experiments 
probing the ionosphere. This was conducted under a S 3 moptic Inter¬ 
national Year of the Quiet Sun (IQSY) program which is directed by 
Dr. Sidney Bowhill. The experiments for which CSL was responsible 
were of a similar nature as the one conducted April 16, 1964, and 
reported in the CSL Progress Report for March, April, and May, 1964. 

The purpose of these propagation experiments is to measure the elec¬ 
tron collision frequency profile of the D region for three different 
conditions; i.e., before, during, and after sunrise. 

All three propagation experiments were successfully conducted 
in the closed loop mode, in which a servomechanism controls the extra¬ 
ordinary wave transmitted power, as explained in the previous progress 
report. In addition to all the signals which were telemetered from the 
rocket to the telemetry stations in the April 16 experiment, the 
receiver automatic gain control voltage was included in the three 
July 15 experiments. The availability of this information, together 
with rocket antenna and receiver calibration data, allows a determina¬ 
tion of ordinary wave electric field intensity versus height to be 
made. It also allows corrections of the differential absorption pro¬ 


file for nonlinear receiver characteristics. 
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Since the sun's ultraviolet light is the main ionizing agent 
for the D region, conducting a series of rocket probes at precisely 
known intervals before and after sunrise provides for the measurement 
of the photoionization buildup due to the sun. Hence, one should observe 
very little differential absorption and Faraday rotation before sunrise 
as compared with after sunrise. The actual differential absorption and 
Faraday rotation data, versus time after launch, for rocket firing 
14.144, during which the earth's shadow was up at about 330 Km, and for 
rocket firing 14.146, at which time there was no earth shadow at the 
launch area, are shown in Figures 1.1 and 1.2, respectively. The dif¬ 
ferential absorption data for flight 14.146 does show a higher trend 
than that for flight 14.144, but not very clearly. The reason for the 
large fluctuations is not presently known. Fluctuations of this magni¬ 
tude were not observed in the differential absorption data for flight 
14.143 fired April 16, 1964, The Faraday rotation data, on the other 
hand, shows clearly the expected dissimilarities between the two flights. 
In flight 14.144, the total Faraday rotation was approximately three 
cycles, whereas in flight 14.146 it was about 13-1/2 cycles. A puzzling 
fact is that the height at which the extraordinary wave was reflected 
in flight 14.143, as indicated by the abrupt rise in the corresponding 
differential absorption and Faraday rotation curves (82 seconds after 
launch, which corresponds to an altitude of about 50 miles) appears 
lower than that for flight 14.146, whereas the opposite should prevail. 

Detailed data analysis is being continued on all four 
rocket firings in order to extract the maximum data. By means of 






Time after launch (seconds) 
Figure 1,1. Absorption Data, 





Faraday rotation (degrees) 
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Figure 1.2. Faraday Rotation Data 
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statistical analysis, Faraday rotation data is being extracted with a 
resolution of about one degree. In this way, a slight decrease with 
height in Faraday rotation before its larger monotonic increase may be 
concluded from the data shown in Figure 1.3 for flight 14.143. This 
data has been corrected for phase shifts which were introduced by the 
electronic circuitry used to extract the data from the telemetry tapes. 
Another circuit is being designed which should eliminate these phase 
shifts. In Figure 1.3 there is a very strong correlation with the 
rocket's precession about its trajectory path, as shown by the sketch 
of magnetic sensor envelope variation versus time. 

A second van installation containing all the ground-based 
equipment except the antenna is being completed for use aboard an 
aircraft carrier. One of the firings aboard the carrier may take 
place at the earth's magnetic equator, at which the rocket trajectory 
would be transverse to the earth's magnetic field. For this transverse 
case a set of mutually perpendicular linearly polarized waves, instead 
of oppositely circularly polarized waves, would exhibit different 
indices of refraction with the resultant differential absorption and 
Faraday rotation. The use of mutually perpendicular linearly polarized 
transmitted waves places a greater burden on the rocket receiving 
antenna, which must be circularly polarized in order to receive both 
linearly polarized waves regardless of azimuth. Circuits to do this 
are now being tested. 


D. Skaperdas 
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Figure 1.3. Faraday Rotation 
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2. SURFACE AND ATOMIC PHYSICS 

Surface: Atomic: 

F. Propst H. Frauenfelder 

R. Retherford F. Franz 

F. Steinrisser 
L. Simonelli 
T. Cooper 

G. Tibbetts 
T. Piper 

2.1 _ Introduction 

In the previous progress report, an experiment to study the 
angular distribution of electrons ejected by ions and by electrons 
striking a solid surface and an experiment to study the adsorption- 
desorption kinetics of gases at solid surfaces were described. The 
design of these two experiments is essentially complete and the apparatus 
is being constructed. It is possible that both of these systems will 
be assembled during the next quarter. 

An instrument designed to study the scattering of electrons 
from solid surfaces has been designed. It is felt that this technique 
is potentially very powerful in that a large number of processes can 
possibly be studied by it. In addition, it may also be possible to 
study ion emission by electron impact using this system. Since high 
resolution (0.1 V) is possible, these experiments may prove to give sub¬ 
stantial data amenable to theoretical analysis—a sorely needed commodity 
in the area of surface science. This system is described in greater 
detail below. 

The apparatus used in the study of Auger electrons (these 
experiments and results have been presented in several previous progress 




reports) has been modified, and work in this area is continuing. In 
addition, the scattering of ions from solid surfaces will be investi¬ 
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gated in this system. 

2.2 _ Rotary-Motion Feedthrough 

The ultrahigh vacuum, rotary-motion feedthrough^ described in 

the previous progress report has been constructed. The output shaft 

2 -10 

was lubricated with dry M 0 S 2 powder. A vacuum of about 5x10 Torr 

was obtained using a small oil pumped system. The exterior bearings 

were water cooled during bakeout. Pressure increments from 2x10 Torr 
“8 

to 5x10 Torr were observed corresponding to angular speeds from 1 RPM 
to about 50 RPM. 

It is expected that the vacuum can be improved by removing 
the exterior bearings and using a clamp for bakeout instead of water, 
cooling. This will allow the bellows and output bearings to be baked 
at a higher temperature. 

T. Cooper 

2.3 Electron Scattering at Solid Surfaces 

2.3.1 Introduction. Indirect evidence that electrons lose 
energy in discrete amounts in interacting with the surface region of 

3 

solids is given by studies of the Auger electron energy distribution 
(caused by slow He ions incident on solid surface) for surfaces with 

^Progress Report for March, April, May, 1964. Section 2.2, 
page 26; Angular Distribution of Auger Electrons . 

2 

The M 0 S 2 powder was obtained from the Moly-Kote Company. 

^F, M. Propst and E, Lilscher, Phys. Rev. 132 , 1037 (1963), 




9 


varying adsorbate coverage. Studies of the photo-electric electron 

energy distribution emerging from surfaces of varying adsorbate coverage 

4 

yield similar evidence. An energy loss mechanism has been used to 
partially explain the energy distributions obtained. However, little 
of positive nature can be learned of this, and other loss mechanisms 
which may be involved, by studying the energy distribution of electrons 
emerging from the surface of a solid which is stimulated by either the 
Auger or photo-electric process. The electron energy distribution 
which is produced within the solid by either the Auger or the photo¬ 
electric process is very broad, the width of the electron energy dis¬ 
tribution being of equal magnitude to the maximum electron energy 
involved. Thus, any energy loss mechanism which may be present will 
not l» resolved and can only be detected by the shifts which they cause 
in this broad energy distribution. 

The electron energy loss mechanism invoked to explain the 

difference between the Auger distribution for a clean W surface and 

3 

that of a W surface with a gas adsorbed is that, when electrons origi¬ 
nating from the Auger process reach the surface of the solid, some 
stimulate electrons in the ground state levels of the adsorbed gas to 
jump to higher levels. In this process, the energy of the Auger elec¬ 
tron is lowered by a discrete amount U, the energy which is taken up by 
the electron associated with the adsorbed gas atom or molecule. The 
net visible effect of this upon the Auger electron energy distribution 
of a clean surface is a decrease in the distribution at the high energy 

4 

C. N, Berglund, Tech. Rep. No. 5205-1, Center for Materials 


Research. 
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end and an increase in the distribution in the lower energy region. 
Further, since all Auger electrons which excite electrons of the ad¬ 
sorbed gas must be displaced in the energy distribution to a point of 
energy at least the amount U below their original energy, the distribu¬ 
tions for clean and gas-adsorbed surfaces should be of the same shape 
(one should be able to get one from the other by simply using the cor¬ 
rect multiplying factor) from the high energy cut-off, Eq, down to 
energy Eq-U. This is found to be true experimentally. Also, it is 
indicated that the multiplier necessary to make the adsorbed-gas Auger 
distribution coincide with the clean-surface Auger distribution in the 
range from E^-U to E^ is directly proportional to the adsorbate coverage. 

The same mechanism can be used to explain the difference 
between the photo-electric distribution for a clean and an absorbate 
covered surface. 

There are many energy loss mechanisms other than the elec¬ 
tronic transitions of an adsorbed gas which also make discrete electron 
energy loss transitions possible. However, most of these cannot be 
"turned on" and "turned off" by adsorbing and flashing-off an adsorbable 
gas as can the electronic transitions of an adsorbed gas, so that their 
shifting effect on the broad energy distribution of the Auger and photo¬ 
electric process cannot be detected', nevertheless, it is quite likely 
that they are present. Some of the energy-loss mechanisms which could 
cause electron energy losses are the following: 
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A. Studies by Corner^ on surface ionization of adsorbed CO 
indicate that molecular vibration levels can be excited. Since these 
are separated by only roughly 1/10 volt, they will be even more diffi¬ 
cult to "see" than the adsorbate electronic levels which have level 
separations of roughly 10 volts. Suppose that the adsorbate was stimu¬ 
lated by a "monochromatic" electron beam of total half width 1/45 volt 
incident on target surface and that the resultant "reflected" beam was 
energy analyzed with the same 1/45 volt resolution, then one hopefully 
should be able to see a set of peaks separated from the incident energy 
by the vibrational energy level separations (since the over-all energy 
resolution is less than 1/25 volt which is less than 1/10 volt). 

B. In semiconductors, studies of surface conduction have 
established the existence of "potential troughs" at the surface of the 
semiconductor. The existence of the surface states in these potential 
troughs, their width, and their exact location could be further estab¬ 
lished by using the same technique suggested above for the study of 
the vibration levels of adsorbed CO. 

C. Acceptor and donor levels provide energy levels which 
make discrete electron-energy loss transitions possible. Even though 
these levels are in the bulk material, this could possibly be studied 
as suggested in A and B. 

^R. Comer and L. Swanson, J. Chem. Phys. 39 _, 2833 (1963). 
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D. The first, second, and etc., filled bands below the con¬ 
duction band of a metal provide energy levels which can cause electrons 
in an incident beam to lose energy. Several types of electron energy 
transition occur if the incident beam has enough energy to excite elec¬ 
trons from these filled bands: (a) the electron initially excited may 
be set free or end-up in the conduction band; (b) an electron from the 
conduction band will fall into the hole exciting an Auger electron from 
the conduction band; (c) etc. 

2.3.2 Experimental Apparatus. The experimental apparatus 
proposed and now under construction to study these energy loss mechan¬ 
isms is a set of electrostatic energy analyzers and the necessary re¬ 
lated apparatus. The electrostatic energy analyzers shown in Figure 2.1 
consist in essence of portions of concentric cylinders with a differ¬ 
ence of potential between. The magnitude of the potential between the 
two portions of cylinders determines the energy of the electrons that 
can pass through the analyzer. The width of the slits at either end 
of an analyzer determines to a large degree the energy resolution of 
the analyzer. The optimum energy resolution of each of the analyzers 
being constructed is 2AV/V = 1/45. 

A single crystal of tungsten will be the first target used. 

A careful energy analysis of the electrons "reflected" from the surface 
of the target as it is bombarded with "monochromatic" electrons for 
both clean and gas adsorbed targets should allow one to make some direct 
and meaningful statement concerning the energy loss mechanisms at work. 





13 


Electrostatic 



Xtal Target | | 

Electron Source Current Detector 


Figure 2.1, Schematic of Electron Scattering Apparatus. 
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2.4 _ Total Cross-Sections for Scattering and for Ionization of 

Barium Atoms by Low Energy Electrons 

Very little data exists in the literature on the properties 
of barium as a plasma constituent. No reference has been found, for 
instance, to measurements of collision cross-sections or to excitation 
functions. This lack of information about barium and other alkaline 
earth elements may be ascribed to their low vapor pressures as compared 
to the alkalies and to zinc, cadmium, and mercury, all of which have 
been studied in some detail. However, advances in materials technology 
have made it possible to contain barium vapor, and an experimental 
barium-plasma device was reported recently by Kennedy, Shefsiek, and 
Talaat. This is a thermionic energy converter in which barium ions 
are employed to neutralize the electronic space-charge. The results 
show that the ionizability of barium vapor, although not as great as 
for cesium vapor, is sufficient to provide the ions needed for neutrali¬ 
zation in a thermionic converter. The degree of ionization is sub¬ 
stantially greater than predicted from the ionization potential and 
the Saha-Langmuir formula. This suggests that the metastable states 
in barium, for which spectroscopic evidence exists, must play an 
important role in producing the ions. However, this role will remain 
essentially speculative until certain fundamental data are obtained. 

The immediate objective of this investigation is to obtain experimental 

”a. J, Kennedy, P, K. Shefsiek, and M, E, Talaat (Martin 
Company, Nuclear Division),"Research and Development for Barium Vapor- 
Filled Thermionic Energy Conversion Technology," Summary Report 
MND-2933-2 (May, 1964), 
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information on total collision cross-sections for scattering and for 
ionization of barium atoms in collisions with low energy electrons. 

Plans have been made for extending the investigation to include excita¬ 
tion cross-sections of the metastable states at an appropriate time. 

An atomic beam method will be used in the study. The essence 
of the method is illustrated schematically in Figure 2.2. A narrow 
beam of barium atoms is produced by an oven and collinating slit. The 
beam is cross-fired by a stream of slow electrons of known energy, 
thereby producing scattering of both electrons and barium atoms, and in 
addition excitation and ionization when threshold energies are exceeded. 

The total cross-section for ionization is the simplest from 
the viewpoint of detection. A suitable positive ion collector can be 
used to count the ions. 

The total cross-section for scattering by slow electrons will 
be determined by a method similar to that used by Rubin, Perel, Bederson 
and Englander in measurements of cross-sections of alkali atoms,^ If 
the acceptance width of the beam detector is less than the beam width, 
scattering collisions deplete the detected beam intensity, thereby per¬ 
mitting the collisions to be counted. A surface ionization detector 

g 

described by Hay will be used to detect the barium atoms. 

^Phys. Rev. 117 , pp. 151-158 (1960), and Phys. Rev. 128 , 
pp. 1148-1154 (1962). 

^R. H, Hay, "The Nuclear Magnetic Moments of Ba^^^ and 

Ba^^'," Phys. Rev. pp. 75-86 (1941). 
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Close control of surface conditions is essential for the func¬ 
tioning of the surface ionization detector and for the detection scheme 
under consideration for metastable barium atoms. The investigation of 
surface effects in these detectors is an important aspect of the program. 
Consequently, the study will be carried out under ultrahigh vacuum 
conditions. 

Preliminary plans have been made for the entire study. 
Construction of components of the apparatus will be started shortly. 

R, Retherford 
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3. COMPUTER RESEARCH APPLICATIONS 

R. M. Brown 

S. Fenves 

3.1 _ Introduction 

This group is concerned with the design and applications of 
digital computers for information processing, primarily in the areas 
of real-time and non-numerical operations. The present work concerns 
the development of an experimental time-sharing system and studies in 
the applications of computers in experimental physics. 

3.2 _ Logical Design 

The designer of digital data processing equipment invariably 
finds himself dealing with Boolean expressions which describe the func¬ 
tions and operations performed by the equipment. More often than not, 
these expressions may be simplified, yielding simpler and more economi¬ 
cal circuits. Although there are many well known techniques for simpli¬ 
fying Boolean expressions, even the best of these techniques can be 
quite time consuming, especially if a large number of variables is 
involved. 

To eliminate this time-consuming operation, a program has 
been written for the CSX-1 which can be used to simplify Boolean expres¬ 
sions containing up to eight variables. The program uses the iterated 
consensus technique first described by Shannon and Mills. 


R. Trogdon R. Jenks 

J. Stifle 
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A description of this technique may be found in "Theory and 
Design of Digital Machines," by Bartee, Lebow, and Reed, Chapter 4, 
Sections 4.1 to 4.8 (McGraw-Hill, 1962). 

J, Stifle 

3.3 _ The DIOG Experimental Time-Sharing System 

During the previous quarter, hardware has been installed in 
the CSX-1 to handle the communication between 7094 and the CSX-1 for 
the proposed DIOG system (see Progress Report for March, April, May, 
1964), This equipment employs five CSX-1 input-output channels; three 
channels provide for data transfer to and from the 36-bit-word memory 
of the 7094 while the other two channels are used for control informa¬ 
tion, This equipment makes the DIOG system compatible with the existing 
7094 system written by the Digital Computer Laboratory, eliminating the 
need for any special program interface between the 7094 system program 
and the input-output channels. 

The Phase Zero DIOG Program has been completed and checked 
out by simulation. In Phase Zero, the user operates in one of two 
modes. In the executive mode, a dialogue is carried out between the 
user at the console flexowriter and the CSX-1 for the purpose of speci¬ 
fying and activating a desired functional program in the 7094. Four 
directives or commands may be given to the DIOG control program in the 
executive mode: "LOGIN" initiates a dialogue; "CALL" enables the user 
to specify an on-line functional program; "TERMINATE" is used to termi¬ 
nate the dialogue with the previously called functional program; 

"LOGOUT" closes out a dialogue. 
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Following the establishment of a program, the user and CSX-1 
enter the program mode where the CSX-1 transmits messages back and forth 
between 7094 and the user. Here messages are accumulated by the CSX-1 
and monitored for one of two signals: (1) a "carriage return" charac¬ 
ter which causes the transmission of the previously entered line to 
the 7094; (2) a "STOP" character flexowriter which produces a transfer 
of the dialogue back to the executive mode. Further details of the 
DIOG system are given by internal memo D2. 

Initial Phase Zero checkout with the 7094 will involve only 
data transmissions between the two computers with the monitor system 
in the 7094 simulated in the CSX-1. Further Phase Zero checkout awaits 
changes in the existing 7094 executive system which will allow the 7094 
to operate either in the on-line mode or the conventional batch 
processing mode. 


R. M. Brown 
R. Jenks 

3.4 _ SMP Bubble Chamber Data Processing 

The SMP measuring console has been moved to the Physics 
Building and a communication line to the CSX-1 computer installed. 

In the coming months three more units will be added and multiplexed 
to the CSX-1. The measuring program, IMSCAMP, is in the process of 
modification to permit measurements on film from the Argonne 30" chamber 
concurrent with measurements on Berkeley 72" chamber film. Significant 
improvements have been made in the accuracy of the data processing and 
the incorporation of console diagnostic procedures. 


R. M. Brown 
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4. SYSTEMS 

Control Systems 
J. B. Cruz, Jr. 

S. D. Agashe 
D. Gieseking 

T. E. Mueller 
T. Murata 

W. R. Perkins 
R. A. Rohrer 
D. Snyder 
G. Tahim 
R. Werner 


Switch. Systems 
S. Seshu 
K. E. Batcher 

A. Dervisoglu 

G. Dodd 

B. Elliott 
W. Mayeda 
K. Onaga 
A. Paul 

J. A. Resh 

H, Y, Chang 


4.1 _ Generation of Trees without Duplication 

One of the important problems in the implementation of computer 
methods of network analysis is the generation of the trees of a linear 
graph without duplication. Further, in order to be able to compute the 
sign of a tree (which has to be done in the analysis of active networks) 
the trees must be generated by replacement of one branch at a time. A 
tree generation scheme satisfying these requirements has been obtained 
and has been reported in report number R-220. In this procedure dupli¬ 
cations are avoided by two devices. On the one hand, we arrange the 
branches of the starting tree into an ordered sequence such that any 
leading subset of branches in the sequence constitutes a connected sub¬ 
graph. On the other, we require that any link which replaces a branch 
of the present tree belong not only to the fundamental cut-set of the 
replaced branch with respect to the present tree (which is necessary 
in order that the new set of branches be a tree) but also to the funda¬ 
mental cut-set of the branch with respect to the starting tree. The 
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replacement is also ordered left to right in the ordered sequency repre¬ 
sentation of the starting tree. It has been shown that this scheme 
generates all the trees of the graph and no tree is generated more than 
once. 

The detailed procedure and the proofs of the assertions above 
are to be found in the cited report. 


W. Mayeda 
S. Seshu 

4.2 _ Quasi-Linear Sequential Machines 

The investigation on the realizability of quasi-linear 
sequential machines has been completed. A necessary and sufficient 
condition was obtained for a given connection matrix (or a state dia¬ 
gram), with the given encoding of states, to be realizable as a quasi- 
linear sequential machine making use of minimal number of feedback 
lines. Also obtained was a canonical realization which enables one to 
realize any arbitrary state diagram as a quasi-linear sequential machine 
using a sufficiently large number of memory elements. These results 
are discussed in detail in Report R-216. 

The invariant properties of the state diagram of quasi-linear 
sequential machines, with respect to the change of coding of states, 
are being studied. 


4.3 


H. Y. Chang 


Self-Diagnosis 

Work has been started on a problem of self-diagnosis; the 


object of study is the laboratory's CSX-1 computer. Initial attempts 
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to devise a bootstrap diagnostic procedure (no external testing device, 
none of the machine assumed error-free £ priori ) have been encouraging. 
The primary tool used is the SEQANALZ (Report R-207) system of programs. 

E. Manning 

4,4 _ Synthesis of Sequential Machines 

A new method for the synthesis of asynchronous sequential 
machines is being developed. The most novel feature of this procedure 
is that it breaks the problem of finding a realization of a given flow 
table into a series of smaller problems in a manner somewhat analagous 
to the synthesis of conventional one terminal-pair electrical networks. 
Further, the procedure is suited to machine computation. 

An algorithm is being developed in connection with this 
synthesis procedure which will attempt to obtain state assignments 
with reduced dependency. Particular emphasis is being placed on obtain¬ 
ing cascade realizations because of the engineering advantages they 
enjoy. 


A. McKellar 


4.5 _ Directed Graph Synthesis 

The study of directed graphs and their application to 
synthesis of switching networks was continued. Some operations on 
a connection matrix which leave the associated output matrix invariant 
were formulated. The relationship between the number of directed 
paths between two distinguished vertices and the nullity of the graph 
is being studied. 


A. Paul 
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4.6 _ Lossy Communication Nets 

A structure of basic saturated cut sets of a lossy communi¬ 
cation net is being investigated. Synthesis procedures of T and T 
matrices of a lossy, bi-complete tree net have been derived. 


K. Onaga 


4.7 _ A-Matrix Synthesis 

Necessary and sufficient conditions for a matrix to be 
realizable as the A-matrix of an RLC network are developed. The RLC 
network is assumed to be non-degenerate or half-degenerate and it is 
assumed to have a connected resistive part. 

The resistive part N and the reactive part N of the RLC 

R X 

network are considered separately. It is shown that if there exists a 
realization then the given matrix A can be factored as -AA^ where A is 
a diagonal matrix of positive entries and A^ is a S3anmetric- 
skewsymmetric (hybrid) matrix. The matrix A determines directly the 
values of capacitances and inductances in the network. A technique is 
given by which the terminal matrix of the resistive part N and the 
fundamental circuit matrix of the reactive part N can be obtained from 
the matrix A^. 

It is shown that the given matrix A has a realization with a 
half-degenerate or non-degenerate RLC network which has a connected 
resistive part if and only if the factorization exists and both the 
terminal matrix and the circuit matrix are realizable. It is also shown 
that for a certain factorization if a realization exists, then the 
reactive part is unique within two-isomorphism and the resistive part 
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can always be reduced to a complete polygon in general. The synthesis 
technique applies to the case where inductances and capacitances are 
time dependent but resistances are time independent. These and other 
results are presented in CSL Report R-225, July, 1964. An extension to 
the degenerate case is being studied. 


A. Dervisoglu 


4.8 _ Dasher Synthesis 

The extension of the Dasher synthesis procedure (realization 
of complex transmission zeros by means of cascaded RC networks) to 
zeros in the right half plane in the sector 


tt/ 3 < jarg s| < tt/2 


is now complete. The conditions to be satisfied by the driving point 
function have all been derived as inequality constraints thus simplify¬ 
ing the preparation problem. The section used in the synthesis pro¬ 
cedure is the simplest possible section which is capable of realizing 
transmission zeros in the desired sector. The procedure for converting 
the unrealizable section which results in this process into a realizable 
section using the Baker node insertion process has also been derived. 

The detailed description and justification of the procedure 
have been reported in report R-231. 


B. Elliott 
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4.9 Distributed Network Synthesis 

Most of the time this summer has been spent on the develop¬ 
ment of a digital computer program to implement the time-domain 
distributed-network sjmthesis method originated by Professor Rohrer. 
Although it is not yet completed, the major difficulties appear to 
have been ironed out; of course, evaluation of the iterative scheme 
being employed will have to be deferred until the program is working-- 
I see no path toward an ja priori convergence proof. 


J. Resh 


4.10 _ Coding Theory 

An aspect of coding theory called the strong converse is being 
considered. This result states that if the transmission rate of a 
communication channel is maintained above capacity, the probability of 
error of an optimal code approaches unity as the code word length ap¬ 
proaches infinity (Shannon's fundamental theorem implies only that the 
probability of error cannot approach zero) the strong converse does not 
hold for all channels, as has been found by Wolfowitz. The problem now 
being investigated involves finding the essential feature pf a channel 
which determines whether or not the strong converse should apply. Pre¬ 
liminary investigation has shown that a channel need not be very compli¬ 
cated for the strong converse to fail; for example, a channel whose 
inputs and outputs are binary sequences of length n, with transmission 
errors uniformly distributed between 0 and n, does not obey the strong 
converse. This result is proved, and some estimates of the relationship 





between transmission rate and probability of error are derived, in a 
paper submitted to the IEEE Transactions on Information Theory . 
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R. Ash 


4.11 _ Nonlinear Oscillating Systems 

Various studies of oscillating systems have been completed 
this quarter. Three laboratory reports were issued and a short 
manuscript written during the period. Two of these reports (R-228 and 
R-229) have been submitted to the journals, and the short manuscript 
(see 4.11.4) has already been accepted for publication. Abstracts of 
the work are given below; 

4.11.1 J. K. Aggarwal, "A Study of Nonlinear Second Order 
Systems," R-223. 

Nonlinear second order systems, which are described 
by a pair of coupled first order differential equations, have been 
studied. The behavior of trajectories in the neighborhood of a singu¬ 
lar point and the behavior of trajectories in the large have been 
examined. 


The isolated singularity at the origin for the system 


k = 
f = 


ni nz 

ax + by 

n3 n4 

cx + dy 


where a, b, c, and d are real, and n^^, n^, n^j and n^ are positive inte¬ 
gers, has been investigated. Sufficient conditions for the stability 
and instability of the singular point have been found. 
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The global behavior of the systems 

= -f(x) + g(y) + p(t) 
f = -h(x) 

and 

k = -f(x) + g(y) + p(t) 
y = +h(x) - k(y) + q(t) 

has been investigated. Here f(x), g(y), h(x), and k(y) are polynomials 
of odd degree with leading coefficients positive, and p(t) and q(t) 
are periodic bounded functions of time. Sufficient conditions have been 
found under which the trajectories of the above systems may eventually 
be confined in a subset of (x,y,t)-space, thus giving bounds on periodic 
solutions (if they exist). Further, bounds on the periodic solutions 
have been investigated by finding regions in (x,y,t)-space from which 
all trajectories eventually leave, and into which no trajectories enter. 
Also, other bounding curves for the generalized Rayleight equation have 
been investigated. 

4.11.2 P, J. Ponzo and N. Wax, "On Certain Relaxation 
Oscillations; Confining Regions," R-228. 

Relaxation oscillations described by the generalized 

Lienard equation, x + iJif(x)x + g(x) = 0 (• = d/dt), with [J, » 1, are 

investigated in the phase and Lienard planes. When f(x), g(x), and 

F(x) = f(u)du are subject to certain restrictions, a number of 
•^o 

analytic curves can be obtained in these planes which serve as bounds 
on solution trajectories. Piece-wise connection of such bounding curves 
provide explicit annular regions with the property that solution 






29 


trajectories on the boundary of an annulus move to the interior with 

increasing time, t. The Poincare-Bendixson theorem then guarantees at 

least one periodic orbit within such an annulus. It is shown that the 

periodic orbits which are isolated by this means are unique within the 

annulus, hence orbitally stable. The maximum width of the annulus is 
“4/3 

of order [j, ; the amplitude bounds obtained for the periodic solution 

agree favorably with the known amplitude for the specific case of the 

2 

van der Pol equation X + |j,(x -l)x + x = 0. The results are generalized 

to less restrictive f(x), g(x), and F(x) than those first considered. 

4,11.3 P, J, Ponzo and N. Wax, "On Certain Relaxation Oscil¬ 
lations; Asymptotic Solutions," R-229. 

Periodic solutions of the generalized Lienard equa¬ 
tion X + [ji,f(x)x + g(x) = 0 (• 5 d/dt) , with H » 1, are investigated 

in the phase and Lie^nard planes. Certain comparison equations are ob¬ 
tained by modifying the functions f(x), g(x), F(x) = f(u)du, so that 

o 

the resulting equations may be integrated to within an error of order 

-2 

|j, . The comparison solutions are used to approximate the solution tra 

jectories of the original equations, and, in particular, the periodic 

orbits. The result is an analytic description of the trajectories, in 

-2 

both planes, to order ij, , The as3miptotic form of the amplitude and 

-2 

period, for a periodic orbit, is obtained with errors of order |j, and 

|j, respectively. For the particular case of the van der Pol equation 
2 

S + |j, (x -l)x + X = 0, the amplitude agrees with the expressions pre¬ 
viously obtained by other workers, while the period does not. The re¬ 
sults of a numerical study of the expression for the period given here 
support the present work. 
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4.11.4 P. J. Ponzo and N. Wax, "On the Periodic Solution of 
the van der Pol Equation." 

The general results of R-229 have been particularized 
to the van der Pol equation. Our analytical results do not agree with 
other formulas for the period. We have computed the period numerically 
using the laboratory's CDC 1604 and established the validity of our 
results. 


J. Aggarwal 
P. Ponzo 
N, Wax 

4.12 _ Finite-Time Stability and Maximum Principle for Systems 

Satisfying a Generalized Lipschitz Condition 

In CSL Report R-215, differential controlled processes 
satisfying a generalized Lipschitz condition were considered. Such 
processes are a generalization of time-varying linear processes. Some 
conditions regarding the existence of a solution of the differential 
equations describing the process and existence of varied solutions, and 
upper bounds for variations of the solution were pointed out. 

It was shown that if a system satisfies a generalized 
Lipschitz condition in the state variables, it is finite-time stable 
with respect to the initial state. If it satisfies a generalized 
Lipschitz condition in the control as well, it is finite-time stable 
with respect to the control also. 

A new derivation of the Pontryagin Maximum Principle and the 
Transversality Condition as necessary conditions for a local extremum 
(and, hence, for a global optimum) was given. This follows a calculus 
of variations approach but uses the Lagrange multipliers in a different 
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way. It was shown that if a generalized Lipschitz condition in the 
state variables and the control is satisfied, and also another condi¬ 
tion is satisfied, for sufficiently small but arbitrary variations of 
the optimum control, the terminal conditions of the problem can still 

t 

be met. 


S. Agashe 


4.13 _ Linear Time Lag Systems 

The concept of complete controllability has been extended 
to the case of linear time invariant systems with time lags. Through 
this the study of the relationship between the frequency domain 
analysis of time lag systems and the optimal control theory is being 
made. The sensitivity of the multivariable time lag systems has been 
investigated from the point of view of Cruz and Perkins.^ 


G. Tahim 


4.14 Stochastic Optimal Control 

Initial investigations, reported previously, have been 
extended to a class of continuous-time, multivariable, noisy control 
systems. 

Studies have focused on a particular type of functional, and 
have led to the definition of related performance criteria. The nature 
of the corresponding optimal control has been investigated numerically 
for specific examples, and analytically for a sub-class of these problems. 


J, B, Cruz and W. R, Perkins, "A New Approach to the Sensi¬ 
tivity Problem in Multivariable Feedback System Design." Presented at 
J.A.C.C., Stanford (June, 1964). 
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Analytical studies are being broadened, numerical techniques 
are being formulated, and results are being generalized to discrete¬ 
time systems. 


M. Sain 


4.15 _ Optimum Time Lag Systems 

2 

It has been shown by N, N, Krasovskii that the optimum 
feedback controller for a time lag system described by 


x(t) = A(t)x(t) + B(t)x(t-h) + C(t)u(t) 

and minimizing 

T 

J = I* [x^(t)Qx(t) + ^^(t)Ru(t)] dt 
•^o 

has the form 

U = Q'(t)x(t) + P(t,T)x(t+T) dT 

•^-h 

where Q'(t) and P(t,T) are matrices. 

The numerical solution for a and 3 is difficult. A gradient 
technique is being developed to find a and 3 by forcing a closed loop 
solution rather than the standard open loop solution used in conven¬ 
tional gradient techniques. 


T. Mueller 


Krasovskii, N, N,, "Analytic Construction of an Optimal 
Regulator in a System with Time Lags," PMM Vol. 26, No. 1 (1962). 
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4.16 _ Optimum Control with a Prescribed Structure 

The following two problems are being investigated. Given 
the optimum control problem 


X = ^(x,u,t) 

final T given 
J(u) = J h(x,u,t) dt 


x(t ) = X , 


(a) find a function ^(x,u,t) [at least approximately] such that 


u = ^(x,u,t), u° to be chosen 


yields an optimum control; 

(b) find matrices A(t) and B(t) [at least approximately] such that 


u = A(t)x + B(t)u 


yields an optimum control. 

Translation of the book "Fondamenti di Calcolo delle 
Variazioni" (Fundamentals of the Calculus of Variations) by L. Tonelli 
is under way. This book introduces certain concepts that have not yet 
been exploited well and that make the Calculus of Variations a proper 
sub-area of functional analysis. 


S. Agashe 




4.17 


Pulse-Width Modulation Systems 

The problem of optimal open-loop control for a system con- 


34 


taining a pulse-width modulator has been solved by the method of the 
steepest-ascent method, and possible utilization of other methods are 
presently under investigation, 

D. Kirk 


4.18 _ A Finite Z-Transform Pair 

A finite z-transform pair has been developed. This pair may 
be thought of as the discrete analog of the finite Fourier transform 
(Fourier series). Consider a periodic sequency of period N, i.e., 
f(n) = f(n+N), Then the transform pair is given by 


f(n) = Y \ 

k=“oo 

N-1 

I f(n)exp(-jnk ^), 
n=0 

Applications of this transform to the parameter variation 
problem in discrete systems and to the steady-state analysis of sampled 
data systems with periodic inputs have been made. 


J. B, Cruz 
W. Perkins 

4.19 _ Optimal Rendezvous Strategies 

A study of optimal rendezvous strategies for cooperative 
controlled systems was completed. The term "rendezvous" is used to 


mean that the systems approach each other so that the velocities and 
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accelerations are the same when the positions become the same. A solu¬ 
tion to the single-input problem, using the calculus of variations, is 
given when the performance index is the integrated-square-error plus 
weighted control. By rewriting the problem in terms of the error coor¬ 
dinates, the rendezvous conditions are replaced by the stability condi¬ 
tions on the Euler-Poisson equations. The stable trajectories are 
found by factoring the coefficient matrix of the transformed Euler- 
Poisson equations. By using the stable trajectories, it is possible 
to express the optimal control laws in terms of the states of the 
systems. The factorization required to obtain the stable trajectories 
is, in general, much easier than trying to satisfy the transversality 
conditions and then deriving the control law. 

A solution to the multi-input problem using Pontryagin's 
Maximum Principle is given for the same performance index. The per¬ 
formance index is rewritten in terms of the stacked state vectors. The 
rendezvous conditions are adjoined as constraints to the performance 
index using Lagrange multipliers. Details of this work are found in 
Report R-218. 

D. Gieseking 


4.20 _ Graph Theory 

The work on unistor graphs started in the last quarter has 
been completed. 

A unistor is an oriented edge connecting two vertices in a 
graph. A flow exists in the oriented edge with a value equal to the 
product of the edge weight and the weight of the initial vertex. The 
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unistor is characterized by flow, potential, admittance, and 
directedness--four properties also found in linear equations and 
electrical networks. 

It was shown that a graph of unistors can represent a matrix 
and the determinant and cofactors can be evaluated by using directed 
trees and directed 2-trees. The signs of the determinant and cofactors 
are calculated during the expansion--a separate sign calculation rou¬ 
tine is not necessary as it is in some other methods. An extension of 
the directed tree evaluation method gave a technique for obtaining 
solutions to a set of simultaneous linear equations. 

A technique was established for representing in a node¬ 
admittance form an active network by a unistor graph. This graph was 
then'used to determine the values of the unknown voltages of the 
original circuit. 

A vertex potential for a communication net was defined and 
used to obtain upper and lower bounds on the maximum flow between two 
vertices of the net. A special class of nets was found for which the 
terminal capacity can be calculated using unistor graphs. 

A complete description of this work is found in Report R-217. 


G. Dodd 
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5.1 _ Introduction 

The purpose of the PLATO project has been to develop an 
automatic computer-controlled teaching system of sufficient flexibility 
to permit experimental evaluation of a large variety of ideas in auto¬ 
matic instruction including simultaneous tutoring of a large number of 
students in a variety of subjects. The PLATO system differs from most 
teaching systems in that the power of a large digital computer is avail¬ 
able to teach each student since one such computer controls all student 
stations. The project work has fallen into three categories, no two of 
which are wholly separate from each other: (1) development of the tools 
for research; (2) learning and teaching research; (3) provision of a 
prototype for multi-student teaching machines. 

5.2 _ PLATO III System Equipment (PLATO Hardware) 

During this quarter work continued in the development and 
construction of circuitry required for the realization of a 20 student 
station teaching system. 

Student station circuitry constructed to date includes all 
circuitry required for the operation of two student stations. A con¬ 
tinuing exemplary effort by various technical sections of the laboratory 


*CSL Consultant 
**Summer Assistant 




38 


dealing with circuit packaging and fabrication suggests that at least 
a total of 10 student stations will be operable some time in November. 

Present interface circuitry has been expanded to include 
the CSL CSX-1 computer as an alternative computer facility to the CDC 
1604. Checkout of the interface was facilitated by a compact general 
engineering routine written during this quarter by Mr. G. Frampton. 

Development continues on special circuitry which will update 
present circuitry or provide special system facilities. Included in 
this circuitry is transistor deflection, power control, master keyset, 
and master video switch circuitry. 

B. Voth 

5.3 _ PLATO III Computer Programming (PLATO Software) 

5.3.1 The Resident Program for CATO (CATORES). CATORES has 
been functional for some time as an input-output monitor for machine 
language teaching programs communicating with the PLATO III equipment. 
During this quarter, a careful examination and re-evaluation of the 
program has taken place with respect to efficient use of the PLATO III 
equipment, with respect to features suggested by the users as desirable, 
and with respect to ultimate compatibility with CATO-compiled programs. 
As a result, the program was almost completely rewritten, with many 
additions and extensions. 

Three different methods of recording student data on magnetic 
tape are now available. For code-cheeking, one can now output student 
data to the printer or to the typewriter as each key is pushed. 
Elaborate checks have been established to avoid time-consuming storage 
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tube operations for a student if they are unnecessary. The plot and 
selective erase routines have been expanded and generalized for easy 
use with CATO-compiled programs, and several subroutines have been 
added expressly for use by CATO programs. In addition, the addresses 
of the subroutine entries and the structures of several internal lists 
have been modified to meet restrictions imposed by CATO. 

Finally, three manually-entered routines have been added. 
PROGSAVE allows an entire teaching program to be written on an auto- 
loadable magnetic tape. SPECTRE allows a student's actions to be played 
back in real time, half time, or double time by executing data read 
from the student data tape. RESTART allows a program to be continued 
from the same point at which it was stopped (in case a lesson runs more 
than one day, or in case of equipment failure). The program uses the 
data on the data tape to reproduce the exact situation present in 
memory and on all student screens at the time stoppage occurred. 

A comprehensive machine language engineering program was 
also written in conjunction with the revision of CATORES. This program 
was used extensively to check and to evaluate the changes to CATORES 
as they were made and is now used to facilitate the maintenance of the 
PLATO III equipment. 

A, Hanson 
M. Walker 

5.3.2 DOPEREA, PLATO "dope" is defined as the record of a 
student's operation of a PLATO keyset. The number of each key pushed 
is stored on magnetic tape along with its associated time, mode, and 
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student number. The stored records may be analyzed at a later date in 
any appropriate fashion. A machine language FORTRAN subroutine, 

DOPEREA, was written which reads and translates the PLATO "dope" tapes 
to binary numbers fast enough to allow smooth reading. 

M. W, Walker 

5.3.3 The CATO (Compiler for Automatic Operations) System. 
CATO (Compiler for Automatic Teaching Operations) was completed during 
the quarter. The system contains three major portions; 

FORTBIN - An expanded and modified version of the FORTRAN-60 
algebraic compiler for the 1604. 

LOGICOMP - A logical compiler which constructs a teaching 
logic interpretable by CATORES from a vocabulary of 10 directives. 

SYSTEMS TAPE - a set of modified subroutines and programs to 
allow easy use and modification of the system and all programs compiled 
with it. 

S. Singer 

A. Hanson 

5.4 _ Plasma Discharge Display Tube Research 

The purpose of the research on the plasma discharge display 
tube is to develop a less expensive replacement for the present PLATO 
storage tube system. The feasibility of initiating the discharges from 
outside the tube was tried with excellent success. Controlling the 
discharge from outside the gas environment not only simplifies construc¬ 
tion by allowing the electrodes to be in the external environment, but 
also effectively adds series impedance to each bulb which helps isolate 
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each cell from its neighbors. The freedom from adjacency initiated 
firing will be investigated next quarter, since the firing of adjacent 
cells should no longer be a problem. 

The cell which was used to demonstrate the feasibility of 
externally controlled discharges was constructed in the following 
manner. Three glass cover plates 25 mm square by 0.15 mm thick were 
prepared for use in a vacuum system. On two of these a gold strip 
25 mm by 0.127 mm by 100 X was deposited on the side that was to be 
on the outside of the cell. The third slide had a 0.127 mm hole 
drilled in its center. A 0.05 mm thick mylar spacer to channel gas 
into the cell was placed on either side of the plate with the hole 
(cell) and the other plates on either side with the gold contacts out¬ 
wards, The whole assembly was connected to glass tubing with epoxy 
and filled with neon. 

R. Willson 

5.5 _ PLATO Learning and Teaching Research 

5.5.1 Inquiry Training (REPLAB). A four week workshop was 
conducted in June in the College of Education of the University of 
Illinois to instruct 28 elementary teachers in the use of scientific 
inquiry classroom techniques. A demonstration class of 22 pupils drawn 
from the local community was held in conjunction with the workshops. 

The pupils used the PLATO program, REPLAB (described in an earlier 
progress report), as one of the lessons in inquiry. The lesson is 
based on a bimetal strip experiment. No statistical analysis was made 
of the work of the summer students. For comparison with classroom 
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behavior many of the workshop teachers monitored the students while 
they were using PLATO, watched reruns of student performance, or 
received copies of the "dope" sheets (student records). Interested 
teachers were also permitted to try the lesson themselves. The rerun 
of student performance was made possible by a computer program, 

SPECTRE, which simulates students using PLATO by reconstruction of the 
operation from the "dope" record stored on magnetic tape during a 
PLATO lesson. 

R, Suchman 

E. R. Lyman 

5.5.2 Teaching A B C*s. Four three-year-old children spent 
two short sessions on PLATO using ALPHABAT, a program designed to 
teach the letters of the alphabet. The PLATO equipment was modified 
to include a semi-automatic audio system. The students were first 
shown a letter on the screen, told its name, and asked to match it. 

Later on, they were not shown the letter, but were simply asked to 
press the key named. The children's reactions and progress gave in¬ 
sight into methods of effective programming for, possible simplifications 
of, and additions to the PLATO system for use in primary and pre¬ 
primary education. 

A. Alpert 

5.5.3 TEXT TESTER. During the past three months, a new 
program called TEXT TESTER has been developed. This program is de¬ 
signed to be used in testing new textbooks. The program provides for 
reproduction of text materials on slides with the insertion of student 
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answers from the keyboard as in PLATO tutorial logics. The student is 
permitted wide freedom in the use of the text; therefore, he is allowed 
to turn pages and to answer questions at random. He is permitted but 
not required to check answers using a judger which compares students' 
answers to lists of stored correct answers. Provision is made for the 
insertion by students of paragraph-form comments about the text mate¬ 
rials. Provisions for testing the student at various points throughout 
the text are also included. 

A major feature of this program is the inclusion of "author 
modes" by which the author may insert or modify correct answers, and 
insert comments which are intended to help students at specific points 
in the program, or which are intended to elicit comments from the 
students. 

TEXT TESTER is the first program ready for use with the 
PLATO Compiler. Code-checking of the TEXT TESTER program will be 
carried on in September followed by trials of the program with students 
later in the fall. 

J, Easley, Jr. 

H. Gelder 
W. Golden 
B. Wilson 

5.5.4 PROOF. During this quarter, work has continued on the 

preliminary programming and advanced planning of the program for the study 

of mathematical problem solving (first described in CSL Report R-185), 

J. Easley, Jr. 

H. Gelder 
W. Golden 
B. Wilson 
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6. VACUUM INSTRUMENTATION 

D. Alpert L. Simonelli 

W. Schuemann F. Steinrisser 

6.1 _ Mass Spectrometer Studies 

As mentioned in the last progress report, we have obtained a 
working Davis and Vanderslice mass spectrometer. Its only problem was 
that the gain of the photomultiplier was about 100 times smaller than 
expected from factory specifications. The reason was probably an 
accident during an earlier bakeout, when a glass tube broke, and the 
system was exposed to air at 250°C. We tried to increase the gain of 
the photomultiplier by introducing a few grains of CsNO^ into the system 
and baking it up to 420°C, as many people had done with success. All 
we observed was a reduced sensitivity of the instrument which can only 
be explained by distortion of the glass envelope causing a misalignment 
of the mass spectrometer. 

We replaced this mass spectrometer tube by a new one, and 
baked the zeolite trap and the rest of the system separately. A valve 
between the system and trap was closed when baking the trap, and the 
latter was immersed in liquid nitrogen when baking the rest of the 
system with the valve open. This time, the instrument met the factory 
specifications. 

A Bayard-Alpert gauge with one carbon-free filament and a 
molybdenum grid purer than usual was added to the system. With the 
valve to the diffusion pump closed, and the Bayard-Alpert gauge running 
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at 10 mA, we observed an equilibrium pressure of 4x 10 Torr, the 
residual gas being mainly He which diffuses through the glass walls, 
expecially the hot wall of the Bayard-Alpert gauge. Cooling this wall 
with an air blower reduces the equilibrium pressure to about half the 
value given above. 

We observed a reduced CO partial pressure in our system com¬ 
pared to the one with a Bayard-Alpert gauge with regular electrode 
material. More experiments are necessary before we can definitely say 
that this is caused by the absence of carbon in the filament and/or 
grid. 

The softness of this purer electrode material gives mechanical 
troubles: the filament is bent towards the grid and, without special 

precautions, touches it easily, and the molybdenum grid is strongly de¬ 
formed during outgassing. 

When the Bayard-Alpert gauge was outgassed with the valve to 
the pump closed, we do not find the cracking pattern of the pump oil 
in our mass spectrum. This can be explained by a complete trapping of 
the backstreaming oil by the zeolite trap. We also find an unexplained 
dependence of total and partial pressures on the outgassing current. 
Further experiments are necessary before any conclusions can be drawn. 
6.2 _ Photocurrent Suppressor Gauge 

Three modifications of the photocurrent suppressor gauge were 
tested this quarter. These three plus the previous model are shown for 
comparison. P„o is identical to the original suppressor gauge except 

JO 

that the size of the suppressor had been increased somewhat to give im¬ 
proved suppression characteristics. See Figure 6,1. 





Figure 6.1. Modifications of Photocurrent Suppressor 
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The general purpose of ^40 continue to improve 

the suppression characteristics and simultaneously to decrease the 
size of the aperture through which the x-rays pass, before striking 
the collector, in order to reduce the unsuppressed x-ray limit. A 
number of interesting conclusions were reached as the result of the 
tests of these two tubes: 

1. The suppressor voltage which is necessary to suppress 
the x-ray photocurrent can be reduced from 300 volts for the original 
gauges at 1x10 Torr to 20 volts by the improved geometr}^ of and 

^ 40 * 

2. The x-ray effect due to photoelectrons leaving the 

suppressor and going to the collector, as the result of x-rays which 

are reflected from the collector and then strike the suppressor, was 

increased significantly for tubes P^g and P 4 Q. In the original sup- 

-13 

pressor gauges this x-ray limit is in the low 10 Torr region while 

-12 

in P^g and P^^ it is in the low 10 Torr region. This increase is 
undesirable. 

3. The sensitivity of the original suppressor gauge was 

20 Torr ^ with a 3/4" aperture in the shield. P^g had a 3/8" aperture 
and a sensitivity of 15 Torr P^q had a 1/4" aperture and a sensi¬ 
tivity of 12 Torr This indeed was surprising and indicates that 
the focusing in this gauge is extremely good. 

With these results, P^^^ was designed. The following advan¬ 
tages of this tube were hoped for; 
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1. A significant reduction in the unsuppressed x-ray limit 
due to the geometry changes. Two 3/8" apertures 1-1/8 inches apart 
are to be used for a shield rather than a short tube. This gives a 
reduction in the x-ray flux coming through the shield by about a factor 
of 10 to 100. The two apertures for the shield rather than a tube 
eliminates ai^ x-rays which would be reflected into the collector region 
due to reflection from the tube, thereby reducing the unsuppressed 
x-ray limit. 

i 

2. A new suppressor geometry (along with the changes in the 
shield) makes it impossible for any x-rays to hit the suppressor 
after one bounce in any place where the resulting photoelectron would 
go to the collector. This should place the x-ray limit, due to photo¬ 
electrons going from the suppressor to the collector, at least a 

4 6 

factor of at 10 to 10 below the unsuppressed x-ray limit, depending 
on the geometry. 

3. The new geometry is much easier to fabricate and align 
than the previous suppressor gauges. 

4. The envelope of is fabricated of aluminosilicate 
glass to make possible lower pressures. 

Preliminary data on this tube indicate that the unsuppressed 

-12 

x-ray limit is in the low 10 Torr region, an improvement by a 
factor of 100. The sensitivity, however, is down by a factor of 5. 
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The preliminary data indicate that a study of the focusing 
properties of these tubes is desirable for two reasons. 

1. To make the design of the suppressor gauge more efficient. 

2. To examine the possibility of building an unsuppressed 

-12 

gauge with an x-ray limit in the low 10 Torr region. 

This study is now in progress. 

6.3 _ Repairing of Damaged or Unusable Storage Tubes 

The PLATO group at C.S.L. uses an expensive Raytheon storage 
tube in their equipment. For two reasons, these tubes have given the 
group considerable trouble: 

1, The total gas pressure in the tubes is too high, causing 
the formation of an "ion spot" on the screen, sometimes in as short a 
time as 10 seconds. This spot both removes the information on the 
screen in that area and creates a broad white area in the center of 
the student's monitor. 

2. Due to accidents the screen suffers a "burn" which has 
the characteristic of either not storing or not retaining the informa¬ 
tion which is written on it. These "burns" were optically visible on 
the storage screen. 

By using several of the tube elements, and with the help of 
a magnet, it was possible to start a Penning discharge in the tube, 
which pumped the tubes rather well. With this pump operating, a number 
of defective tubes were baked at temperatures ranging from 100°C to 
200°C, sometimes with the filament on, and sometimes off. All tests 
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successfully lowered the equilibrium pressure during normal operation 
to a sufficiently low pressure that the ion-spot formation time is of 
the order of 3-5 minutes. This processing successfully eliminates the 
problem of ion spot formation. 

In addition, it was noticed that baking the defective tubes 
seemed to diminish the seriousness of the "burns." Subsequent process¬ 
ing, where the screen was heated to 350°C and the gun to 150°C for an 
extended time with the Penning discharge on, has resulted in the com¬ 
plete repair of one tube, and it has been returned to service. An 
interesting fact is that the optically visible "burn" disappears during 
processing. Further experiments are preceding on the "burn" repair 
problem. 

6.4 _ Impact Ionization Studies 

During this quarter, work has continued on the impact ioni¬ 
zation of adsorbed gases on surfaces due to electron bombardment. 

CO has been put on a moly sample, and its maximum cross sec- 
tion for impact ionization determined to be about 6x10 ions/electron. 
The temperature of desorption for the adsorbed CO is about 800°C, 

When oxygen was placed on a platinum sample, the yield of 

_8 

ions/electron was less than 10 . If the filament of the Bayard-Alpert 

was left on (giving a conversion of O 2 to CO), the yield increased to 
6x10 The gas layer was desorbed at about 200°C. The values for CO 
were 6x10 ^ and 200°C also, indicating that the CO is adsorbed on the 


surface in both cases. 
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The amount of CO on the surface seems to be a small fraction 
of a monolayer because (1) flashing the sample gives a very small in¬ 
crease in the pressure, and (2) the time required to saturate the effect 
when gas is being introduced is much less than that required by a mono- 
layer. Further studies using an all-metal system are in progress. 
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7. PLASMA PHYSICS 

C. Mendel 
H. G. Slottow 
M. A. Smith 
R. Hosken 


W. Carr 
J, Crowder 
B. Hicks 
T. Lie 


7.1 _ The Boltzmann Equation for a Shock Wave 

We developed an algebraic theory of shock wave structure in 
the previous quarter. A draft of a report on the algebraic theory was 
finished and sent to Professor Yen at AVCO, who will be one co-author, 
for his comments. The theory, within the limitations imposed by the 
nature of the basic assumptions, covers all Mach numbers and any 
monatomic gas. 

There have been two applications. The first application was 
to the calculation of shock structures for the full range of Mach 
numbers, from one to infinity, for a gas of elastic spheres. Direct 
comparisons between the results of the algebraic theory, the Mott- 
Smith theory, and Navier-Stokes theory are possible. A second appli¬ 
cation of the algebraic theory was in the analysis of production runs 
with the Monte Carlo program. This use of the theory indicates that 
the Monte Carlo calculations are converging too slowly, and that one 
high-order boundary condition is not being satisfied with adequate 
accuracy. 

In an attempt to study convergence of the Monte Carlo program, 
without interference from the statistical fluctuations, we ran fifty 
iterations with a fixed collision sample at each station. These runs 
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were made with 10 collisions/sample, five stations, and a Mach number 
of 2.5. No clear convergence was indicated, possibly because there is 
an interaction between the rate of convergence and the large statisti¬ 
cal fluctuations corresponding to these small samples. A similar but 
less extensive run for a larger sample also did not clearly show ade¬ 
quate convergence properties. In this more accurate calculation, how¬ 
ever, the period of the long-term fluctuations during the iterations 
was decreased by a factor of two or more. 

In view of these several indications of convergence troubles, 
and the recommendations of Dr. Mott-Smith regarding methods of inte¬ 
gration, it seemed appropriate to redesign the integration method, A 
new method has therefore been derived which imposes a least-square fit 
of the derivative of a pol 3 momial, for each station, to the Monte Carlo 
derivatives df/dn at each interior point of the shock, and requires 
exact matching of the polynomial to the boundary values of f. Proper 
choice of the degree of the polynomials will force a certain smoothing 
of the "noisy" values of the derivatives calculated by Monte Carlo. 
Fitting in the manner described to the interior values of the derivative 
and to the end values of the function makes step-wise integration, with 
its danger of instability, unnecessary. We believe that a second step 
toward improving the integration method will make fitting of a polyno¬ 
mial for each velocity bin subject to auxiliary conditions expressing 
the fact that three moments of the derivatives vanish at each station 


in the shock wave. 
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While these revisions of the program are being made, we shall 
undertake a very simple calculation for a "pseudo-shock" that was sug¬ 
gested by Arnold Nordsieck as a test of his Monte Carlo method. 

B. Hicks 

7.2 _ Oscillations in Non-Uniform Plasmas 

In the previous quarter we reported observations of incoherent 
scattering of microwaves by a high frequency mode of a finite non-uniform 
plasma. Although a number of papers have appeared during the last two 
years that discuss certain features of these modes, the whole subject 
is still in its infancy.^ ^ 

In the following, we outline a method that permits us to 
treat one special case of a non-uniform plasma in considerable detail, 
and that provides a suitable basis for the investigation of more gene¬ 
ral cases. 


We consider a one-dimensional plasma enclosed between two 
infinite, specularly reflecting walls at x = 0 and x = L. The behavior 
of the plasma is described by the linearized Vlassov and Poisson 
equations 


^^1 5^1 e 

^ ir - m ^1 


3f 

n(x) ^ 


= 0 , 


3E 

-— = -4TTe Ff^ dv . 
Sx J 1 


( 1 ) 

( 2 ) 


The term E Sfi/Bv can be shown to be small, and it will be neglected 
o i 

henceforth. As boundary conditions, we require 


^The papers cited here are listed in the Bibliography (§7.3), 
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f (0,v) = f (0 ,-v) , 
f(L,v) = f(L,-v), 

E(0) = E(L) = 0. 


(3) 

(4) 


Following Weissglas we introduce the functions 
= f(v) + f(-v) and F = f(v) - f(-v). The boundary condition for 
the distribution function now simply requires F = 0 at the boundaries, 
The functions f and E are assumed to have a time dependence of 
Equations (1) and (2) then become 


2 ^2 - ^ M 

F + ^ ^ V ' + 2‘i —■* n(x) E. = 

2^2 mto lov 

(JJ Bx 


0 , 


jV V dv}. 0. 

o 


If no external driving field is present, equation (6) implies 

CO 

4nei 


"l = 


(B 


J F V dv , 


(5) 


( 6 ) 


(7) 


For F~ and Ej^, we now choose the following expansions which satisfy the 
boundary conditions; 


F = E Fj^ sin(Trkx/L); E^^ = E Ej^ sin(TTkx/L), 


( 8 ) 


For n(x) we choose the special profile 


n(x) = + '^cos2trx/L). 


( 9 ) 


It will be shown later that this choice is less restrictive than it 
might appear. Inserting (8) and (9) into (5) and (7) results in the 
following recurrence relation for Ej^; 
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( 10 ) 


with 


u) L +00 Of /3v)dv 
_ 1 . P (• o _ 

®k “ TTk J ^ (JU - TTkv/L 


For k = 1 we have to observe that E . These are two 

infinite sets of homogeneous equations, for odd and even indices. 

For a solution to exist, the determinants must vanish. This 

requires 


0 = 


and 


€ 2 ^ " 2 2 ." 
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1(63-1) 
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|( 63 - 1 ) 


( 11 ) 


0 = 
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0 




|(C6-1) 


( 12 ) 
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In order to evaluate these determinants we divide the 


row deteminant then becomes 


0 = 


26 , 


v(e^-l) 


- 1 


2e, 


^(63-1) 


2e, 


v(e^-l) 


• • • • • • 9 i 


and a corresponding expression holds for the even set. The determinant 
is of the form 


a^ 1 0 0 


D = 


1 1 


0 1 1 


I • • • • 


We divide the first row by and subtract it from the second row. This 
yields 

I 

1 /a. 


D = a. 


1 

0 

0 


a^-l/a^ 


0 0 

1 0 


• • 
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Proceeding in this fashion we obtain 

D ■ “l<Vl^aj)(a3 - —^)(a^--)... . 

=3 • a^-l/aj 

By calculating the partial numerators and denominators of these con¬ 
tinued fractions it can be shown that all terms cancel except for the 
numerator of the last continued fraction. Hence 


D = lim A , 
n 

n-«o 


where A is the numerator of 
n 


a 

n 



If this continued fraction vanishes, we may write 


0 = 



(13) 


This is the dispersion relation in the non-uniform case. 

A particularly simple expression for the continued fraction 
arises if we replace by its value derived from the moment equations 
for a warm plasma 


O) 




- (TTkv^/L)^ 


(14) 


Inserting this expression into (10), we notice that the recurrence 
relation for the is identical with that for the coefficients of the 


Mathieu functions se (x,q) , if we put 

n 
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q 



D 


nvT 

U) 

P 


Hence, we conclude that E obeys a Mathieu equation 


E" + (a - 2q cos2ttx/L)E = 0 . 


The eigenmodes of oscillation are then given by the corresponding 
values of a and q, belonging to the Mathieu function se^(x,q). The 
corresponding eigenfunction E^ is the Mathieu function se^(x,q). 

The qualitative behavior of these eigenvalues is sketched in 
Figure 7.1. We notice the asymmetry for positive and negative values 
of q corresponding to the two different density profiles (Figure 7.2), 

The eigenvalues can of course be calculated exactly by using 
the correct expression for the dielectric constant and solving for the 
roots of (13). Computer calculations in this direction are in progress. 
Likewise, the eigenfunctions corresponding to the modes can be obtained 
numerically. 

The knowledge of the eigenfunction and eigenvalues provides 
a starting point for the calculation of more difficult density profiles 
by perturbation theory, because the near degeneracy of the eigenvalues 
for V = 0 has been removed, 

E. Jackson 
M. Raether 
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Figure 7.2. Plasma Density Profiles 
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8 . HIGH MAGNETIC FIELD SUPERCONDUCTORS 

C. B. Satterthwaite R. Ries 

J. 0. Kopplin R. Klingbiel 

W. 0. Gentry M. G. Craford 

R, Bernstein 

8.1 _ Introduction 

This group is primarily interested in studies which contribute 
to an understanding of the persistence of superconductivity to very 
high magnetic fields. Studies include materials preparation, investi¬ 
gations of the properties of Type II superconductors, and the properties 
of superconducting specimens of very small dimensions. 

8.2 Superconductivity in Thallium-Antimony Alloys 

In an attempt to find a Type II superconductor with a transi¬ 
tion temperature near the lambda-transition of liquid He, the supercon¬ 
ducting properties of a series of thallium-antimony alloys have been 
studied. The series contained alloys of composition between 0 and 30 
atomic percent of antimony. 

The samples were prepared in an inert atmosphere and heat 
treated for several hours to attain homogeneity. The transition tem¬ 
peratures were obtained from magnetization studies using the vibrating 
sample magnetometer. From these measurements were also obtained the 
critical fields He, Hcj^, and Hc 2 characterizing the Type II supercon¬ 
ductivity in these alloys. Measurements of the residual resistance 


ratio were also made. 
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Residual resistivity measurements indicate that in the neigh¬ 
borhood of 13.5 atomic percent Sb an ordered structure must occur. 

Also in this region of composition there is a distinct rise in the 
transition temperature, though the maximum occurs between 15 and 18.5 
atomic percent. 

The alloys were Type II superconductors, but they were far 
from ideal in that their magnetization was highly irreversible, and 
there was appreciable trapping of flux. 

This program has been terminated for the time being with the 
completion of a Master's thesis, the termination of Mr. Gentry's tenure 
as an industrial research fellow, and his return to the Fansteel 
Metallurgical Corporation. A report will be forthcoming. 


W. Gentry 


8,3 _ Superconducting Tunneling 

The anomalous current peaks in tunneling characteristics of 
Nb-Nb oxide-In junctions, reported in the previous QPR, have been 
studied in some detail. This structure has now been observed in a 
number of junctions composed of a variety of superconductors, and it 
is now believed to be characteristic of junctions exhibiting a DC 
supercurrent resulting from either a tunneling supercurrent (the dc 
Josephson effect) or a low critical current superconducting short. 

In situations where two superconductors are separated by a 
thin insulating barrier, we know from Josephson and others that one 
can transfer Cooper pairs at zero-voltage from one side to the other 
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without change of energy. At finite voltages, however, the transfer 
involves a change of energy 2eAV per pair, where AV is the voltage dif¬ 
ference between the two sides, and the transfer is allowed only if 
there is an absorption mechanism for photons of this energy. 

A current at a finite voltage in excess of the quasi-particle 
tunneling current is then evidence of such an absorption process and 
is a manifestation of what has come to be known as the AC Josephson 
effect. 

In Figure 8.1 the peak structure in the current-voltage curves 
is shown relative to the normal two-superconductor tunneling 
characteristic. 

To resolve the structure, which represents a whole series of 
negative resistance regions, we resorted to the circuit shown in 
Figure 8.2. With the parts of the circuit indicated immersed in 
liquid He, it was possible to provide a very low impedance source for 
the tunneling current. The load line, which is numerically equal to 
Rl + R 2 , could be made as low as a few milliohms, thus allowing us to 
resolve very sharp current peaks. 

Tunneling current peaks at finite voltages have been observed 
under a number of circumstances, and it is possible to make some 
generalizations about them. 

1., The current peaks at finite voltages always accompanied 
a super-current peak at zero voltage and were of comparable magnitude. 
Similar structure was observed whether the zero-voltage current was a 








Figure 8.2. Circuit for Measuring Tunneling Current. 
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Josephson tunneling current or a low critical-current superconducting 
short. This is not unreasonable, since the two are, in many ways, the 
same. It is clear that if a voltage can be sustained across either, 
the transport of pairs must be accompanied by an absorption process. 

The shorts could be easily distinguished: (a) they had no magnetic 
field dependence, (b) the currents at zero-voltage were greater than 
the theoretical maximum for the DC Josephson effect, and (c) there ap¬ 
peared to be an appreciable normal conductance parallel to the tunneling 
current. 

2. Variation of temperature affected the structure in two 
ways. First, the number and amplitude of the peaks increased with 
decreasing temperature. This corresponds to the enhancement of the 
Josephson effect or the increase in the critical current of a short 
with decreasing temperature. The second effect is that the voltages 
at which the peaks occur are shifted uniformly with variation of tem¬ 
perature. The shifts in voltage with temperature variation for both 
Nb-NbO -In and Nb-NbO -Sn junctions, relative to the transition tem- 
peratures of In and Sn, are shown in Figure 8.3. 

3. In cases where the current peaks accompanied superconducting 
shorts, the structure was independent of magnetic field, but when they 
accompanied Josephson tunneling the amplitude of all peaks were periodic 
in magnetic field and with approximately the period of the zero-voltage 
peak. They were, however, not necessarily in phase with it or with 

each other. This is evident from Figure 8.4, which shows curves for 
the same junction at a series of magnetic fields. ' 









Current (Arbitrary units) 


















71 


In most of the junctions there appeared to be series of 
peaks at multiples of a fundamental frequency. For example, the 
series of peaks indicated with circles in Figure 8.4 is spaced at 
multiples of approximately 10.7 microvolts; the first three peaks 
occur quite accurately at this spacing. Those peaks marked with 
triangles form another series at 20 and 40 microvolts. 

It appears that the structure we observed is very similar to 
the steps in the tunneling curves observed by Fiske. These have been 
analyzed by Fiske and his co-workers as standing wave modes in the 

junction itself. Treating the junction as a transmission line which 

has a capacitance per unit length, determined by the thickness and di¬ 
electric constant of the insulating barrier, and an inductance per 

unit length, determined by the penetration depths, they have been able 
to calculate standing wave frequencies of the magnitude observed. One 
gets from a treatment of this kind a temperature dependence of the 
voltage at which the peaks should occur, resulting from the temperature 
dependence of the penetration depth. 

The voltage at which the peaks should occur is inversely 
proportional to where is the penetration depth of Nb 

and X is the penetration depth of the other superconductor. In or Sn 
in the present case. Therefore, the ratio of the voltage at any tem¬ 
perature t to the limiting value V^> as the reduced temperature t ap¬ 
proaches zero, is V/V^ = vt 

Inserting values for the penetration depths of Nb and In and their 
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temperature dependences one obtains the theoretical curve shown in 
Figure 8.3. The theoretical curve for Nb-Sn is not appreciably dif¬ 
ferent from the one shown. 

M. Craford 
R. Peacock 
R. Ries 

C. Satterthwaite 

8.4 _ Micro-micro Voltmeter 

A trial version of a cryogenic, low level DC voltmeter has 
been constructed and is being tested. Sensitivity and other operating 
characteristics are being studied. 

The voltage sensitivity of this instrument is primarily due 
to its unique input circuitry. The input of the voltmeter consists of 
a small solenoid which is vibrated inside of a set of pickup coils. A 
sensitive AC amplifier and lock-in detection system, connected to the 
pickup coils, detect the presence of any DC current in the solenoid. 

The entire input circuit is operated in a liquid helium bath, allowing 
the use of superconducting wire in the solenoid and pickup coil. 

The input impedance is therefore purely inductive, so that 
extremely small DC voltages can be detected. The voltage must be large 
enough compared to the input inductance so that the measurement can be 
made in a reasonable length of time, and the source impedance must be 
below some minimum level, however. 

Preliminary measurements indicate the sensitivity to be 
quite high, although considerable difficulty is encountered from 
spurious signals. The main spurious signals are apparently caused by 
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poor shielding of both the earth's field and the AC fields at the instru¬ 
ment's carrier frequency. One microampere in an inductance of 200 jiH 
is presently detectable, although superimposed upon spurious signals 

equivalent to a 100 |j.A input. With this sensitivity, and source resist- 
-7 -12 

ances of 10 ohm or less, one picovolt (10 volt) would be detectable 
in about three minutes, and would be well above the noise in twenty 
minutes, 


8.5 Crystallization of Nb^Sn 


R. Ries 


Considerable progress has been made toward growing Nb^Sn 
single crystals of sufficient size for the study of their properties. 
Crystals with flat faces up to 2 mm on an edge have been obtained. 

A total of ten runs has been completed using the Nb^Sn 
crystal-growing furnace. The duration of these runs varied from 19 
hours to 160 hours. In all cases a niobium-tipped rod was inserted 
into the Nb-Sn solution to provide the site for nucleation of the crys¬ 
tals. This produced polycrystalline starts as numerous crystals were 
obtained from each run. Figure 8.5 shows the niobium tip and Nb^Sn 
crystals as taken from the furnace following run number 5. The solu¬ 
tion temperature, in this case, was held at 1286°C and the initial tip 
temperature was 1220°C. The solution was stirred magnetically and the 
tip was in contact with the solution for a total of 43 hours. 

Figure 8.6 shows a single crystal after being isolated and removed 
from the tip. 




Figure 8. 


Figure 8. 


. Nb^Sn Crystals as Withdrawn from Furnace, 
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Correlation between crystal structure and size with the 
various furnace parameters has been limited because of failure and/or 
necessary modification of certain furnace parts. All crystals have 
been grown from a molten Nb-Sn solution held at a temperature of ap¬ 
proximately 1290°C with a maximum temperature variation of 2°C during 
the crystal growing period. Withdrawal or raising of the niobium tip 
during the crystal growing period at any rate approaching .0075 inches 
per hour was detrimental in that only a high number of very imperfectly 
structured crystals were obtained. Rotation of the tip at 12 rpm 
resulted in the growth of larger crystals with fewer structural 
imperfections. 

A schematic diagram of the furnace is shown in Figure 8.7. 

The Nb-Sn solution is contained in a high purity alumina crucible 
supported in the hot zone of the furnace by a lava pedestal. A pt/pt 
137o Rh thermocouple adjacent to the silicon-carbide heating element 
was used to control the furnace temperature, and thermocouples inside 
the niobium tip and in a ceramic tube inserted into the solution were 
used to measure the tip and solution temperatures, respectively. A 
slight flow of argon was established and maintained through the furnace 
during the crystal growing period after the furnace chamber had been 
evacuated using a standard vacuum pump. 

The next phase of the investigation will use small Nb^Sn 
crystals for top-seeding the Nb-Sn solution. The problem here is se¬ 
curing the rather small crystals in a manner which will withstand the 
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high furnace temperature. Attention will be given to determining if 
there is a preferred crystallographic orientation for growth of the 
crystals. Also, studies of the purity, structure, electrical and mag¬ 
netic properties of the crystals are planned. 

Work on the chemical analysis of solution samples taken from 

the furnace, for the purpose of determining the solubility of niobium 

in tin in the 1100-1300°C'temperature range, has continued, but final 

results have not been obtained. Considerable time was spent developing 

a cerate titration technique using nickel as a reductor for the stannic 

tin solution, only to cast doubt on the suitability of the nickel 

reductor for the quantitative determination of tin in this system. 

J. Kopplin 
R. Klingbiel 
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9. HIGH VOLTAGE BREAKDOWN 

D. Alpert T. Casale D. Lee 

E. M. Lyman H. Tomaschke W. Edelstein 

9.1 _ Critical Field for Tungsten Electrodes 

The data for broad-area parallel-tungsten electrodes in 
ultrahigh vacuum have been presented in the form of Fowler-Nordheim 
plots. The slope of these plots yields the enhancement factor P by 
which the average or gross electric field V/d at the cathode must be 
multiplied in order to obtain the local field present at the tips of 
the whiskers responsible for the field emission current. The final 
result for the average value of the critical field Ej^ at breakdown is 

6.5 + 1 X 10^ V/cm, which agrees very well with the results of Dyke 

1 2 3 4 

et al. , Boyle, Kisliuk, and Germer, Gofman et al. , and Tomaschke, 

Preliminary results were reported in the CSL Progress Report for June, 

July, and August, 1962, where it was shown that the critical field 

thus obtained is constant within the experimental error for over five 

orders of magnitude of gap spacing. The value of E^^ may thus be 

^Dyke, Troian, Marten, and Barbour, Phys. Rev. 9 ^, 1043 (1953) 
2 

Boyle, Kisliuk, and Germer, Journ. Appl. Phys. 2 ^, 720 (1955) 

3 

Gofman, Protopopov, and Shuppe, Soviet Phys. Solid State 
Physics 2, 1203 (1960). 

4 

H. Tomaschke. PhD. Thesis, University of Illinois, Urbana, 

(1964). 
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regarded as a characteristic property of the electrode material. A 
very important consequence of this result is that the value of break¬ 
down voltage is thus related to, and directly predictable from, the 
observed volt-ampere characteristic of the predischarge field emission 
current when the work function cp is known. Knowing the critical elec¬ 
tric field at which breakdown will occur, one can arrive in a non¬ 
destructive way at the value of voltage at which breakdown will occur. 
Another important result is that under clean conditions, the breakdown 
characteristics for electrodes of a given material are apparently inde¬ 
pendent of the geometry, and of whether the electrodes are single 
crystal or polycrystalline in structure. 

9.2 _ Multiple Points on Broad Electrodes 

In our broad area electrode experiments, the predischarge 
current was found to arise from a number of point emitters on the 
cathode. As the voltage between the electrodes is slowly raised, at a 
current of about one microampere, the electron beams impinging on the 
anode produce sharply defined blue-green pin points of light on the 
anode. The light is very faint and requires dark adaptation to be 

observed. The brightness of the spots increases with increasing cur- 

2 

rent, and they number from 1 to 10 per cm . For an electrode spacing 
of about 0.4 cm (V = 160 kV, I = 150 1 J,A) , their diameter is approxi¬ 
mately 0.01 cm. Although the spot pattern is relatively unvarying over 
a period of hours, certain spots may disappear or flicker and new spots 
appear elsewhere from time to time, either with or without visible 
breakdown, the total number remaining roughly constant, at a given gap 
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spacing. It was observed that a lateral movement of the cathode with 
respect to the anode (maintaining constant gap spacing) caused the 
entire pattern of spots to move across the anode in conjunction with 
the cathode motion. Thus it seems plausible to attribute each indi¬ 
vidual spot to the beam from an individual whisker or protrusion on 
the cathode. The size of the anode spots agrees with that calculated 
by the method of Vibrans^ for the beam spreading from a whisker of the 
dimensions observed by Tomaschke and others.^ The light intensity 
from the spots is roughly that to be expected from the visible tail of 

the bremsstrahlung spectrum of the impinging electrons, as suggested 

7 8 

by Silverman and Vibrans, or from "Lilienfeld" radiation. 

When the high voltage was initially applied to electrodes 
which had not been exposed to high voltage for a period of a day or so, 
fluorescent spots which did not move laterally in conjunction with the 
cathode motion were seen at the anode. Spots were also seen on the 
cathode. These fluorescent spots were of somewhat different color 
than the spots produced by the impinging electron beams. They com¬ 
pletely disappeared in a minute or so, and are presumably due to specks 

^G. E. Vibrans, Lincoln Laboratories Technical Report 308 
(18 April 1963). 

^Little and Whitney, NRL Report No. 5944 (May 20, 1963) and 
Journ. Appl. Phys. 2430 (1963). 

^E. Silverman, Lincoln Laboratory (MIT), High Power Tube 
Program Semiannual Report (December, 1962), 

g 

See, for example; H. Boersch, C. Radeloff, and G. Sauerbrey, 
Phys, Rev. Letters 7, 52 (1961). 
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of fluorescent material adhering to the anode and cathode; their 
disappearance is associated with gas bursts and sharp initial current 
pulses, which also disappear in a short time after applying the high 
voltage. 

9.3 _ Effect of Gas on the Predischarge Current 

Previous experiments were repeated in an attempt to obtain 
a reproducible measurement on the effect of gas in reducing the pre¬ 
discharge current between tungsten electrodes. Although a reduction 
of two orders of magnitude in the current was observed, presumably due 
to the blunting by ion bombardment of the emitter tips, the experiments 
were regarded as unsuccessful, because the results were not reproduc¬ 
ible. The trouble seemed to be due to a masking of the gap emission 
current by an ion current which takes the longest discharge path (from 
flange to flange) in the electrode chamber, and builds up by the 
Townsend mechanism. The apparatus is being modified with guard rings 
to isolate the ion current from the gap current. 

9.4 _ Predischarge Current and Voltage Breakdown between Gold 

Electrodes 

Gold electrodes of purity 99.95 have not given the degree of 
reproducibility in predischarge current and voltage breakdown that were 
obtained with the single crystal tungsten electrodes. Nevertheless, 
straight Fowler-Nordheim plots are obtained, and the enhancement fac¬ 
tors derived from the plots in conjunction with the observed breakdown 
voltages indicate a preliminary characteristic electric field intensity 
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at breakdown, E^, of 6,4+ 1 x 10^ V/cm, which is nearly the same as 
that obtained for tungsten. (0 = 4.8 eV was assumed for goldi) 

9.5 _ Flicker and Rectangular Pulses 

9 

It was previously reported that the field emission current 
from both small and broad area electrodes exhibits "rectangular" fluc¬ 
tuations as if some of the emission sites were bistable in their 
emission properties. The duration of the pulses ranges from micro¬ 
seconds to seconds, or even longer. Then, it was observed that in both 
types of electrodes the light spots on the anode due to impinging elec¬ 
tron beams fluctuate in intensity, or flicker. When observed on a 
modified field emission microscope, the flicker was found to consist 
of the appearance and disappearance of small bright spots within the 
pattern of a single emitter point. 

By focusing the light from a flicker spot onto a photocell, 
it was found that the emission current fluctuates in exact correspond¬ 
ence with the light. 

Several attempts were made to eliminate the flicker. Cooling 
the cathode support wires with a liquid helium bath produced no ob¬ 
servable change in the flicker, nor did cleaning the cathode by heating 
it to temperatures as high as 2400°C. 

In order to obtain a purer material, a tungsten cathode was 
heated in the presence of tungsten carbonyl vapor, causing the tungsten 
carbonyl to decompose and deposit a layer of pure tungsten on the 
cathode. This reduced but did not eliminate the flicker, 

^CSL Progress Report for Dec,, 1963, Jan, and Feb., 1964, p,74. 
Also Progress Report for March, April, May, 1964, p,89. 
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However, with a cathode made of ultrahigh purity tungsten 
(obtained from the Westinghouse Electric Corporation) no flicker what¬ 
ever was observed, even when the voltage was increased to the point of 
breakdown. 

On the basis of this single preliminary observation, one may 
postulate that the flicker is due to the diffusion out of the bulk and 
migration across the surface of the tungsten on the part of impurity 
atoms. As these atoms hop onto and off of an emission site they alter 
the work function and hence the emission current, thereby creating the 
flicker effect. Further work with ultrahigh purity tungsten is in 
progress to check this tentative hypothesis. 
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10. THIN FILMS 


R. N. Peacock 

R. 

E. Zelac 

M. G. Craford 

W, 

P. Bleha 

J. T. Jacobs 

A. 

Jones 


10.1 _ Tin Oxide Films 

The report for the previous quarter discussed transparent 
conducting tin oxide films and gave results of various measurements 
which were made to learn something of the electronic structure of this 
material. Measurements were reported on the temperature coefficient 
of resistance, stability, and the effects of gaseous and vacuum ambi¬ 
ents on conductivity. To learn more about the structure, the transient 
photoconductivity has been measured with visible and with ultraviolet 
light. Analysis of this data is under way at present. Preliminary 
results show that the band gap in this material may be larger than the 
energy corresponding to the shortest ultraviolet light used. 

J. Jacobs 


10.2 _ Hard Superconducting Films 

Attempts made earlier this year to deposit films of vanadium 
and niobium on substrates at room temperature, even under good vacuum 
conditions, did not produce films having normal superconducting proper¬ 
ties. Transition temperatures, when they were found, were far below 
the bulk values. Using heated substrates, films of V and Nb have been 
made recently which have superconducting transition temperatures simi¬ 


lar to bulk material. 
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The method used for heating the substrates, although unusual, 
has given good results. Substrates are customarily heated by placing 
them in contact with a heated metal or graphite block, or by using 
infrared heaters. When these heaters consume large amounts of power, 
they outgas themselves and adjacent portions of the evaporation system, 
and at best do not insure that the substrate has the required constant, 
uniform temperature. The substrates used in making the superconducting 
films were of quartz 1/16" X 3/8" X 1-1/2" with about 2000 % of titanium 
on the back side. This titanium film served as an ohmic heating element 
in perfect thermal contact with the substrate. Outgassing was negli¬ 
gible, and, since the power consumption was small, there was no heating 
of the remainder of the system. A temperature of 600°C was obtained, 
using about 40 V at 0.5 A, or only 20 watts. In testing, temperatures 
as high as 1000°C were possible without destruction of the Ti film. 

The phenomenon of superconductivity has served in this case 
to show the large difference in the structure of films deposited at 
different substrate temperatures. The ability to make good films of 
superconductors will be useful in future experiments. 

W. Bleha 


10.3 _ Epitaxial Films 

The vacuum system to be used for these studies was completed 
at the start of this quarter, and pressures as low as 1 X 10 Torr 
were obtained. Titanium sublimation pumping, when properly used, was 
quite successful in maintaining a low total pressure in the system. 
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In the operation of this system an unusual amount of trouble was 
encountered with copper-gasketed seals leaking after bakeout, and con¬ 
siderable time has been lost for this reason. Other baked systems in 
the laboratory, using the same flanges, have not had this problem. 

Only tin films have been studied so far, and these were 
deposited on cleaved NaCl substrates which had been baked in vacuum 
before cleavage. The aim of these experiments is to produce highly 
oriented single-crystal films. Electron-diffraction and electron- 
microscope examination of some tin films has shown that they consist 
of oriented polycrystalline films at best, and do not have the desired 
uniform structure. It is possible that the use of thicker films, 
higher or lower pressures during deposition, or better outgassing of 
the tin charge could lead to better results. Substrate temperature 
is an important parameter, of course. 

R. Zelac 


10.4 _ Luminescent Films 

The usual luminescent screens in cathode-ray tubes are made 
by a settling process which can be performed only by a specially 
skilled technician, if uniform screens are to be made. At best, they 
have irregularities which limit resolution. Such screens also evolve, 
into the vacuum, various active gases harmful to surface experiments. 
Evaporated luminescent screens, which have been made by a number of 
workers, do not suffer from these difficulties, and would be useful to 
a number of projects in this laboratory. 
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Suitable screens usually cannot be made simply by vapor- 
depositing a normal phosphor. Evaporation may destroy the activator, 
or stoichometry, or both. Heat treatment after deposition can restore 
luminescence. Selection of phosphor materials is important for this 
method. 

Films of CaWO^ were recently deposited and successfully 
activated by heating in air, although the process was not completely 
reproducible. Experiments with this, and with other phosphors, will 
continue. Rare earth oxides, which can have luminescent properties, 
may be promising according to some reports. 

W. Bleha 
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11. COMPUTER OPERATIONS 


H. G, Slottow 
C. Arnold 
S. Chilton 
G. Crawford 
L. Hedges 


J. Knoke 
V. Metze 
R. Nash 
E. Neff 
J. Stifle 


11.1 _ Introduction 

This group is responsible for the development and operation 
of the CSL computing facility. 

11.2 _ CSX-1 Computer 

11.2.1 Operations 

Period; June 3, 1964, to September 3, 1964 
Total running time; 

Ave. per day (7 day week); 

Operational Time; 99% 

Scheduled Maintenance Time; 

Emergency Maintenance Time; 

Breakdown of emergency maintenance time; 

8.0 hrs; replace 25 volt power supply 
0.33 hrs; replace shorted transistor (2N2270) 

0.17 hrs; replace diode that was installed backwards 

J. Stifle 

11.2.2 Modifications 

Eight additional interrupt channels have been added to the 


859.3 hrs. 

9.44 hrs. 
850.8 hrs. 
0.0 hrs. 
8.5 hrs. 


computer. 
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Additional equipment has been added to facilitate communi¬ 
cation between the CSX-1 and the 7094. 

Additional input-output circuitry has been added to allow 
the PLATO experiment access to the CSX-1. 

J. Stifle 

R. Trogdon 

11,3 CPC 1604 Computer 

11.3.1 Operations 

Period: June 1, 1964, to August 31, 1964 
Total running time; 

Ave. per day (7 day week) 

Operational Time; 95.25% 

Preventive Maintenance Time; 4.56% 

Emergency Maintenance Time: 0.19% 

Breakdown of emergency maintenance time: 

2.00 hrs; Intermittent memory failure. 

E. Neff 

11.3.2 Systems Programming 

A new resident system is being prepared for the CDC 1604 
computer. The assembler for the system, now almost complete, contains 
most of the features of FAP. In addition, it contains a more flexible 
expression evaluator, allowing the use of parenthesis, and a compile 
pseudo operation. Descriptions of this and other programs in the 
system will be published as they are completed. 

R, Nash 

S. Seshu 


1049.9 hrs. 

11.41 hrs. 
1000.00 hrs. 
47.9 hrs. 
2.00 hrs. 
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F,0. Box 73 * 

L*xlii(ton 73. IteaaacteMtta 
Atta! Library, A-OK 

I Sylvaola Elaccrlc Fraducta. lac. 
Elaccroaic Syatn* 

Mltbaa LA*. Library 
100 Fint AvaoM 
Ualchn 54. lhaaachaaatt* 




1 Mr. A.A. Lunditm 

Ball Talapbona Laboratoriaa 

Boob 2E-127 
Uhl^ny toad 
Uhippany, Niv Jarcay 

I AFMSC (NDSCO HaJ. F. Mwalcc. Jr.) 
HclloaBn Air Force Baa* 

Hau Mexico 68330 


4455 Genarae Street 
Buffalo 21, dev York 
AitO! J.F. Danond. Librarian 


Sperry Cyroecope Ccopany 

HariM Dtvlaioo Library 

155 Clen Cove toad 


Major UiUiaB Karri* 
BADC (lAtfl) 


1 tne Air Develon*nt Center 
Grlffie* Air Force Baec 

Attn: Docxpaenc* Library 
BAALD 



ColuaOia 


57, Nev York 





tactalle-WFtindk 
Baccelle Haaiorlal Xnaclcuta 
505 Kina Aveme 
Coluobue L, Ohio 

Aeronivcicel Syeceoe Dlvleloa 
NavLaetlon and Guidance Laborecory 
Wilfhc-Faccereon Air Force Baet 
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22314 


KtiMand Bulldtnc 
3045 Coluabla Fikc 
ArllDtCon 4, fIcBinU 


U.S. lieval Ueepone LAoratory 
CcapuCation and Analyeie Laboratory 
Dahlgren, Virginia 
Attn; Hr. Ralph A. Nleoann 


1 Mlno e apolla-Wooeyvall Radiator 
Aeronaiticat Blvtalea 

Mlanapoita 13, Hinaaaata 


Attn: Br. D.F. Elvall 


1 la t p ec tar of Naval Hacarial 

loren of Ship* Tacbaical Napraaaocatlva 
1*02 Itaac Hinnabaha Avaiaw 

20 Activity Supply Ofltcer. 

talldti« 2504, OMrla* Mood Area 

Fat: AccoumAlc Proparty Offlear 
INrkad: For Inat. for E^epioracary Bceearcb 





















